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Abstract 



> 

I/-) ■ The exact expression for Wilson loop averages winding n times on a closed 



contour is obtained in two dimensions for pure U(N) Yang-Mills theory and, 
rather surprisingly, it displays an interesting duality in the exchange n <-► N . 
The large- N limit of our result is consistent with previous computations. 

q; Moreover we discuss the limit of small loop area A, keeping n 2 A fixed, and 

D 
— ■ find it coincides with the zero- instanton approximation. We deduce that 

:^ ■ small loops, both at finite and infinite "volume", are blind to instantons. 

X 

j^ ■ Next we check the non-perturbative result by resumming 't Hooft-CPV and 

Wu-Mandelstam-Leibbrandt (WML)-prescribed perturbative series, the for- 
mer being consistent with the exact result, the latter reproducing the zero- 
instanton contribution. A curious interplay between geometry and algebraic 
invariants is observed. Finally we compute the spectral density of the Wilson 
loop operator, at large N, via its Fourier representation, both for 't Hooft 
and WML: for small area they exhibit a gap and coincide when the theory is 
considered on the sphere S 2 . 
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I. INTRODUCTION 

The non-perturbative structure of non-abelian quantum gauge theories is still a chal- 
lenging topic in spite of a large amount of efforts in this direction. Whereas perturbation 
theory provides a well-established computational tool to describe the weak coupling regime, 
quantitative predictions for the behaviour of the strongly-coupled theory are extremely hard 
to be found. In the last few years much attention has been devoted to derive exact results 
in the supersymmetric case, by exploiting duality properties [|I]|§. Nevertheless we believe a 
"traditional" field theory approach is appealing, most of all in comparing perturbative and 
non-perturbative aspects. 

In particular this approach becomes crucial in the light-front formulation of quantum 
gauge theories 0: though some non-perturbative features are thought to be transparent, a 
consistent framework in the continuum is still lacking. Moreover the relation with the usual 
perturbative equal-time quantization remains unclear. 

Such problems have recently been tackled in the simplified context of two-dimensional 
gauge theories (YM2), taking advantage of the lattice solutions 0. As far as the continuum 
is concerned, in two dimensions the theory seems trivial when quantized in the light-cone 
gauge. As a matter of fact, in the absence of dynamical fermions, it looks indeed free, being 
described by a Lagrangian quadratic in the fields. 

Still topological degrees of freedom occur if the theory is put on a (partially or totally) 
compact manifold, whereas the simpler behaviour enforced by the light-cone gauge condi- 
tion on the plane entails a severe worsening in its infrared structure. These features are 
related aspects of the same basic issue: even in two dimensions the theory contains non triv- 
ial dynamics, as immediately suggested by other gauge choices as well as by perturbative 
calculations of gauge invariant quantities, typically of Wilson loops 0. One can say that, 
in light-cone gauge, dynamics gets hidden in the very singular nature of correlators at large 
distances (IR singularities). 

The first quantity that comes to mind is the two-point correlator. If the theory is 



quantized in the gauge A^ = at equal-times, the free propagator has the following causal 
expression (WML prescription) in two dimensions 

first proposed by T.T. Wu ||. In turn this propagator is nothing but the restriction in two 
dimensions of the expression proposed by S. Mandelstam [0] and G. Leibbrandt [§] in four 
dimensions and derived by means of a canonical quantization in ref. @. In the equal-time 
formulation, the causal behaviour is induced by the propagation of unphysical degrees of 
freedom (probability ghosts), which can be expunged from the "physical" Hilbert space, 
but still contribute in intermediate lines as timelike photons do in the QED Gupta-Bleuler 
quantization. 

In dimensions higher than two, where "physical" degrees of freedom are switched on 
(transverse "gluons"), this causal prescription is the only acceptable one; indeed causality is 
mandatory in order to get correct analyticity properties, which in turn are the basis of any 



consistent renormalization program |L0] . It has been shown in perturbative calculations [1 1 
that agreement with Feynman gauge results can only be obtained if a causal propagator is 
used in light-cone gauge. 

The situation is somewhat different in exactly two dimensions. Here the theory can 
be quantized on the light-front (at equal x + ); with such a choice, no dynamical degrees of 
freedom occur as the non- vanishing component of the vector field does not propagate 

D p ++ {x) = - l -\x-\5{x% (2) 

but rather gives rise to an instantaneous (in x + ) Coulomb-like potential. On the other hand, 
renormalization is no longer a concern for a pure Yang-Mills theory in two dimensions. 

A formulation based essentially on the potential in Eq. (|2]) was originally proposed by 
G. 't Hooft in 1974 |12|, to derive beautiful solutions for the gg-bound state problem under 
the form of rising Regge trajectories. 

When inserted in perturbative Wilson loop calculations, expressions ([!]) and (0) lead 
to completely different results, as first noticed in ref. [f|. The origin of this discrepancy 



was eventually understood in ref. []13|| , where it was shown that genuine non-perturbative 
excitations ( "instantons" ) are necessary in the equal-time formulation (Eq. ([1])) in order 
to obtain the exact result, which in turn is easily recovered in the light-front formulation 
(Eq. (^|)) just by summing the perturbative series. This surprising feature strengthens the 
belief that the light-front vacuum may indeed be much simpler than the one in the equal-time 
formulation, at least in two dimensions. 

In order to gain a deeper insight in the different physics described by either 't Hooft 
or WML prescription, we have decided to study a wider class of loops, i.e. loops winding 
around themselves an arbitrary number of times. Actually, their knowledge appear to be 



intimately related [14|| to the general solution of the Makeenko-Migdal equations at large 



N p5| . Furthermore, in the large- N limit, the spectral density of the eigenvalues of the 
Wilson operator is completely determined in terms of such loops. Remarkably this spectral 
density carries non-trivial information about the master field |16[ of the theory even on the 



plane. Therefore it seems natural to explore winding loops using both the perturbative 
WML prescription and the 't Hooft one in order to understand whether instanton effects 
are still able to distinguish between the two. 

In Sect. II we start by deriving Wilson loops at finite N and arbitrary n on the plane 
from the non-perturbative solution on the sphere S 2 . To our knowledge an explicit formula, 
obtained by decompactification, appears in literature for the first time and extends previous 



results limited to a small number of windings [17,18|. It contains all instanton contributions 



and turns out to correspond to the exact resummation of the perturbative series defined 
via 't Hooft propagator. Once again light-front quantization seems to capture the correct 
vacuum of the theory. On the other hand, by isolating the zero-instanton contribution on 
the sphere, we end up with the perturbative resummation of WML expressed in terms of a 
Laguerre polynomial. Hence the picture is fully consistent with our previous investigations 

mm- 



Having the general expression for a Wilson loop of U(N) with n windings, we are enabled 
to perform different limits covering complementary regimes of the theory. First of all, in 



the large-iV limit with g 2 N fixed, we recover the well-known result of ||T7| , f20| in terms of 
a Laguerre polynomial. Surprisingly enough, taking instead n to infinity and keeping n 2 A 
fixed, the perturbative result of WML is reproduced. We can explain the latter behaviour by 
observing that, having the instantons a finite size, small loops are essentially blind to them. 
This also teaches us, as we will investigate further in the sequel, that a remnant of the weak 
phase on the sphere survives the decompactification limit. In fact one can as well study 
the same limit directly at finite total area and what is found is precisely the zero-instanton 
contribution. The parallelism of the limits N — > oo and n — > oo, reflected by the same 
functional form, has to be ascribed to the symmetry of the general formula for the Wilson 
loop under the exchange N <-> n, up to a rescaling of the area. 

The next step consists in computing the same quantities by means of the perturbative 
expansion. This is not simply a check of the correctness of the previous interpretation, but 
also casts light on the interplay between geometry and colour in the different prescriptions, 
and is done in Sect. III. Whereas for WML the resummation of the full perturbative series 
is a straightforward generalization of the computation performed in [2I| for n = 1, the 



analog with 't Hooft unveils the full non-abelian character of such a prescription. Due to its 
complexity, we have restricted ourselves to the calculation of the Wilson loop with winding 
n to 0(g 4 ): as we will see, at this level already a large amount of classes of diagrams 
need being considered. As a matter of fact, we have realized that adopting WML all the 
diagrams contributing to a given perturbative order entails the same geometrical factor, so 
that only the structure of the colour traces is relevant. Fortunately the sum of the colour 



factor is encoded in the matrix integral of |H]. At variance with WML, 't Hooft prescription 
leads to different geometrical integrals for different classes of graphs, which, together with 
the colour traces, complicate the resummation of the perturbative series. Still, in this 
framework we have been able to describe the mechanism underlying both the limits iV — > oo 
and n — > oo. It is rather amusing to observe how the latter limit precisely selects a class 
of diagram whose geometrical factor is the same as in WML, explaining in this case the 
perturbative resummation with WML. When performing the former limit the situation is 

6 



trivial from an algebraic point of view (only non-crossed diagrams contribute), while the sum 
over geometrical factors is quite involved. We have evaluated it by exploiting the symmetry 
under the exchange N <-> n, which shows a curious interplay between geometry and the 
algebraic invariants. 

Finally in Sect. IV we have applied previously derived results to the computation of 
the spectral function in the large- N limit. In the 't Hooft case the situation is well-known 
fl6| , |22|l , a gap appears when the area of the loop is below a critical value, resembling the 
weak-phase behaviour of the same quantity on the sphere. The novelty comes out with 
WML, for which a gap was expected at any value of the area, since instanton excitations 
should be neglected. On the contrary the gap is present only for sufficiently small loops and, 
in this case, again the density coincides with the same quantity on the sphere in the weak 
phase. 

Conclusions are drawn in Sect. V, whereas technical details are deferred to Appendices 
A and B. 



II. THE WILSON LOOP WITH WINDING NUMBER n 

Our starting point are the well-known expressions || of the exact partition function and of 
a Wilson loop with winding number n for a pure U(N) Yang-Mills theory on a sphere S 2 
with area A 
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d-R(s) being the dimension of the irreducible representation R(S) of U(N); C 2 (R) (C 2 (S)) 



is the quadratic Casimir, A — A, A are the areas singled out by the loop, the integral in 



is over the U(N) group manifold and Xr(S) is the character of the group element U in the 
R (S) representation. We immediately notice that W n = W_ n . 

Eqs. (0), (|4]) can be easily deduced from the solution of Yang-Mills theory on the cylinder, 
using the fact that the hamiltonian evolution is governed by the laplacian on U(N): we call 
Eqs. (|f),((|) the heat-kernel representations of 2(A) and W n (A — A, A), respectively. 

In order to evaluate W n (A— A, A) in the decompactification limit we write them explicitly 
for N > 1 in the form 
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We have described the generic irreducible representation by means of the set of integers 
rrii = (mi, ..., mjv), related to the Young tableaux, in terms of which we get 



c 2 (R) = t-(N 2 -i) + J2 [m- 
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A is the Vandermonde determinant and the integration in Eq. (U) has been performed 
explicitly, using the well-known formula for the characters in terms of the set rrii. 

In this section we will study the partition function (|j) and the Wilson loop ([|) in the 
limit of infinite area, by showing they are dominated by particular representations labelled 
by suitable indices {mi}. Let us suppose N odd, so that the term (N — l)/2 in the Casimir 



operator is integer and can be absorbed in the sum over rrii [] 
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It is now easy to see that the dominant contributions are given by the following set of indices 

N- 1 , 

(9) 



{m*} H 0, ±1, ±2, • • • , ±- 



with all possible permutations, for which the minimum value of the Casimir is reached and 
reads 



c 2 ({m}) 



N(N 2 - 1) 
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(10) 



Finally, by exploiting the symmetry, in the large A limit the partition function becomes 

9 2 A 



Z(A — > oo) = A (mi, ..., tjin) exp 
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We now turn to the Wilson loop (|6|). Thanks to its symmetry, we can always choose k = 1 
and the equation becomes 



W n (A - A, A) = — exp 
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When the previous formula is evaluated in the decompactification limit A — > oo, ^4 fixed, 
that is for {m^} = {m^}, we get 
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After writing explicitly the Vandermonde determinants and taking into account the symme- 
try over the indices {rh 2 , ..., ttin} we are left with a nice formula for W n 



1 If N is even, (JV — l)/2 will be half-integer, but this does not alter our conclusions. 
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Remembering that Wo = 1 and that W n is even in n, in the following we shall explicitly 
consider only positive values of n. 



It is not difficult to show that (|14j) can be conveniently rewritten as 
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that can be combined together to produce the general expression 
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The series is actually a finite sum, stopping atft = n — lorft = iV — 1, depending on the 
smaller one. Some comments are now in order. First of all we notice that when n > 1 the 
simple abelian-like exponentiation is lost. In other words the theory starts to feel its non- 
abelian nature as the appearance of different "string tensions" makes clear: the Coulomb 
law is violated and the combinatorial coefficients in ([L7]) are intimately related, as we will 
see in the following, to the presence of instanton contributions to light-front vacuum. Actu- 
ally, from the sphere point of view, (|i~7|) can be understood as coming from an instantons' 
resummation. Indeed, as first noted by Witten p3 |, it is possible to represent Z(A) and 
W n (A — ^4,^4) as a sum over unstable instantons, where each instanton contribution is as- 



sociated to a finite, but not trivial, perturbative expansion. In ref. ||13||, as confirmed by 



10 



the computation on the torus in |1| , it was shown that if only the zero-instanton sector is 
considered, in the decompactification limit one exactly recovers the sum of the perturbative 
series for n — 1, in which the light-cone gauge Yang-Mills propagator is WML-prescribed 
(pi). Precisely we had for the zero-instanton case 



W{ 0) = — exp 
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Lg (x) being the generalized Laguerre polynomials, reproducing for A — > oo the exact 



resummation of |21 
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This has to be contrasted with the full result coming from (17) 
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giving the expected area-law exponentiation, that can be easily derived directly on the plane 
by resumming the perturbative series in which the light-cone propagator is CP^-prescribed 
according to 't Hooft [Q. Actually, in the large-N limit, ( |T9"D does not exhibit confinement, 
as first noticed in [2T|. As a matter of fact from (JT9]), taking the limit iV — > oo (g 2 = g 2 N), 
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however, this is hardly surprising since W} does not contain any genuine non-perturbative 
contribution, and simply reproduces the Wilson loop behaviour of the weak-coupling phase 
of the sphere ||24j| , where instantons are suppressed p5[ , |2"6| , p2" . 



For a general winding n, Eq. (p!7|) is therefore expected to come out from the resummation of 
't Hooft perturbative series, that corresponds to a light-front quantization of the theory || . 
On the other hand to neglect instantons, and then to send the area of the sphere to infinity, 
is likely to reproduce the WML computation: the perturbative analysis will confirm these 
claims. At the moment we need the zero instanton contribution to (|lTD: it can be easily 
obtained, closely following ref. [T3[], by applying a Poisson resummation to (^|) and taking 
the non-exponentially suppressed terms as g 2 — » 
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For A — > oo it becomes 
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The above result is quite different from ([L7]) and reffects much more dramatically the same 
discrepancy found for n — 1: the string tension is independent from N and there is no 
trace of the non-abelian nature, since different exponential weights do not appear. Actually, 
the behaviour of ~W n (the WML one) is, with respect to the fundamental winding n = 1, 
exactly the same as in the abelian U(l) theory. It is easily proved for U(l) that 

g 2 An 2 ' 



W n (A; 1) = exp 



(24) 



the n windings resulting trivially in replacing the charge g by gn. The WML result satisfies 
the same rule as it is clearly seen from ( ^3|) 



W^(g 2 A;N) = w[ \g 2 n 2 A;N). 



(25) 



In this perspective the exact (i.e. full instantons 't Hooft) result does not seem related to 
any simple-minded reduction U(N) ~ U(1) N , suggested by the abelianization of the theory 
in axial gauges; rather a non-trivial (non-abelian) effect is present there. 

We can further observe another basic difference between Eq. (|IT) and Eq. (p3[): as a 
matter of fact, the former shows an interesting symmetry between N and n. More precisely, 
we have that 



W n {A;N)=W N {A;n) 



(26) 



A - U A 



a relation that is far from being trivial, involving an unexpected interplay between the 
geometrical and the algebraic structure of the theory. Looking at Eq. (|26|) , the abelian-like 
exponentiation for U(N) when n — 1 appears to be connected to the U(l) loop with N 
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windings, the "genuine" triviality of Maxwell theory providing the expected behaviour for 
the string tension. Moreover we notice the intriguing feature that the large- N limit (with n 
fixed) is equivalent to the limit in which an infinite number of windings is considered with 
vanishing rescaled loop area. Alternatively, this rescaling could be thought to affect the 
coupling constant g 2 — > j^g 2 - Eq. (^), of course, does not exhibit such a fancy behaviour; 
nevertheless it is related, in a particular dynamical limit, to the full result, as we will see in 
the following. 

Let us discuss in detail the large- N limit of Eq. (|T7|) : it was obtained long ago by Kazakov 
and Kostov ]17|Q) who solved the Makeenko-Migdal equations in the case at hand. Coming 
back to Eq. flT5|), we introduce the function W n (A; N) from 

g 2 ANn 



W n (A;N) = exp 
that can be expressed through a contour integral 

-. n— 1 



W n (A;N) 



Nn 



k=0 



g 2 A ( . n - 1 
n k 



k\ 



dtF~ n (l-t) 



n—l 



(27) 



(n- 1- k)\2iri J c 

(C is a contour surrounding the origin of the complex plane). The binomial sum can be 
performed and, after reversing the path around the pole at t = 1 and finally changing 
variable (z = 1 + t), we get 



Wn(A;N) 



2-niN 



cxp 



-g 2 A n{n - l) 



in-1 



dz [l-(l + z)ex V (g 2 An/2)f , 

N M • I 28 ) 



ic {1 + z. 

The case n > N is easily obtained from Eq. ([T7|) by exploiting the symmetry between n and 
N (Eq. (ED). 



2 The explicit form for finite N and small n (n = 2,3) can be deduced from |L7J] and [18|. The 
latter used a non-abelian version of Stokes' theorem, and both are consistent with our general 
formula (|T7|). 
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Eq. (|28|) is a nice representation for the Wilson loop with n windings, which is particularly 
suitable to discuss the large- N limit. Actually, rescaling z into -^, we have 



(-l) n 

ZTil 



-g 2 A 



n(n — I) 
4 



dz [N - (N + z) exp (g 2 A n/2] 



n n-l 



C l 1 + f 



2-\N 



yn+1 



(29) 



and then, by taking the limit N — > oo, g 2 = g 2 N fixed, we arrive at 
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from which we can easily recover the Kazakov-Kostov result 
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Using Eq. ( p6f ) we are able to perform another limit, namely n — > cxd with fixed n 2 ^4 
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Eq. ( p2|) exactly coincides with the WML result Eq. (|23| ) : this means that in the small area 
limit (taking n — > cxo in order to have a non-vanishing interaction) we essentially recover the 
zero-instanton, i.e. perturbative, approximation, even from the exact formula Eq. (0). The 
same behaviour will show up when studying the spectral density for Wilson loop eigenvalues 
(see Sect. IV). Once again the reason for the equivalence of Eqs. fl3~2f) and ( |23"1) is founded on 
the non-trivial geometrical aspects encoded in Eq. (|I7D . In fact, we recall that in Eq. ( |23|) 
we have neglected the instanton contributions: we expect that sufficiently small loops are 
unaffected by instantons due to their having a typical length scale measured by their size. 
Eq. ( |3~2"D supports this intuitive argument. 

This conclusion is not limited to the decompactified case: we can go further and prove 
that, even at finite total area A, the exact heat-kernel expression Eq. (|6]) reduces to Eq. ( |22|) 
when n — >• oo Q. Introducing X = n 2 A, we can write Eq. (|6]) as 



3 We need to have n — » oo with n 2 A fixed so that the interaction is finite. 



14 



W n Z = exp 



48 



N{N 2 



exp 



9 2 X 



oo 

Y^ exp 



m,=— oo 



<?M 



AT 



o _ 5^A mi 
2 n 



x A 2 (m 1 ,...,m N )[[( 1 + 

i=2 ^ 



n 



nil ~ m j 



Expanding the product we get 
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We obtain a series in 1/n from the expansion of the exponential term 
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Eq. ( |33"D is truly an expansion in 1/n 2 : in order to understand this point, we simply notice 
that changing rrii — > —rrii we produce an overall factor (—l) k+l weighting the sum over 
m,j. This implies that k + I (and therefore k — I) must be an even integer so as to have a 



non- vanishing result. It is useful to rewrite the sum over m^ as 



1 oo 

nn £ exp 



(fc+i) 



m;=— oo 



9 A 



E 



i=l 



m- 



fe+i 



A 2 (mi,...,m Ar ) Z JJ 



m 



p(i) 



P j=2 



m P (i) - m P(i ) 



where the sum over P is over the elements of the permutation group Sk+x- Let us evaluate 
to the O th -order term: we observe that for I = k 



fe+i 



En 

P j=2 



m 



p(i) 



rap(i) - m P (j) A (mi, ..., m fe+ i) 



^/ P (m lr .,m w ) 



(34) 



the Vandermonde determinant being produced by the common denominator; the quantity 



Y^p/p(mi, ...,m.fc+i) is a polynomial of degree 



fc(fc+i) 



in fc + 1 variables. Since the Vander- 



monde determinant A (mi, ..., m^+i) is totally antisymmetric, the non- vanishing contribution 
to Eq. ( |3"B| ) comes from the totally antisymmetric part of Y^ p /p(mi, ..., rrik+i), implying that 

^/ P (mi,...,m w ) ~ cA(mi,...,m fc+ i). 
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The constant c is easily proven to be 1 by inspection of Eq. fl34|). Hence the contribution of 
order (4,)° to Eq. © is 

"JV-l 



W n Z = exp 



9 2 A a:i ^ 2 
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J2 exp 



m,i=— oo 



iV(^ - 1) 



2 A N 

9 A 



E 

fc=0 



(A r-l)! (-1)* /^A 



(iV-l-Jfe)! k\(k + l)\ 



E 



«!■ 



A 2 (mi,...,mjv) 



Z[A]jjL%Ug*\/2)e-^, 



(35) 



proving that the O th -order contribution is exactly the WML one. It remains to show that 
the divergent contributions (coming from I < k in Eq. (|33"D) are vanishing. To this purpose 
we use again Eq. fl34|) to notice that we can still factorize a Vandermonde determinant in 
the denominator, but now Y^ p fp(mi, ..., rrik+i) is a polynomial of degree less than *■ * ' in 
k + 1 variables. Its totally antisymmetric part vanishes, so that 



(i) 



Jim W n (A -A,A)~- L^ig 2 A/2) exp 



N 



9 2 X 



+ o - 



n- 



This confirms our guess that sufficiently small loops do not see instantons: actually it is 
known that in the strong phase (at large- N) small loops behave as in the weak phase (i.e. 
they do not confine) as n — > oo. Our result goes further, revealing that even at finite iV and 
at finite total area A, instantons can be dynamically suppressed. 

III. THE PERTURBATIVE WML AND 'T HOOFT PRESCRIPTIONS 



We now turn to the resummation of the perturbative expansion for W n . The aim of our 
investigation is not only to check the correctness of the interpretation of the non-perturbative 
results derived so far, but more importantly to observe and explain a curious interplay 
between geometry and colour factors in the two different prescriptions. We will begin with 
WML, for which the resummation of the perturbative series is found through a simple 
generalization of the procedure outlined in [^TJ in the case n — 1. This is due to all diagrams 
entailing the same geometrical factor. We will then switch to 't Hooft prescription and realize 
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the resummation is rather involved since graphs belonging to the same class, according to 
their colour factor, produce different geometrical integrals. Hence, in this respect, light-front 
quantization exhibits its peculiar non-abelian character. 

A. Resummation of the perturbative series defined via WML prescription 



We closely follow the procedure outlined in PH] . In the euclidean space the perturbative 
expansion of W n is 



1 _ _ C P s 2k—1 

W n [A] = l + -J2(-9 2 ) k d Sl x^( Sl )--- " ds 2k x^(s 2k )Ti[G, 1 ... IM2k (x ir --,x 2k )]., 
fe=i ^° ^° 

(36) 

where x fJ, (s), s G [0, 1] parametrizes the contour with n windings. The Lie algebra- valued 
2A;-point Green function G fll ... fl2k (xi, • • • , x 2k ) has to be expressed via the Wick rule in terms 
of the free propagator D(xi — Xj), in the light cone gauge A_ = and endowed with (Wick- 
rotated) WML prescription 

Dt + {x-x') = ^5 abX -±^±. (37) 

Z7T £__ — X_ 

The simplest choice for the contour is a circle wrapping around itself n times. Then it 
happens |21[ that the weighted basic correlator is independent of the loop variables, since 



• / \ • / l\ X +\ S ) X +\ S ) / \2 

X- (s) X- (s — , ; - = 2(mrr) 



X-[s) — x'_(s') 



After that, the integration over the path parameters s±, ■ ■ ■ , s 2k trivially yields l/(2k)\ and 
the task of determining the Wilson loop reduces to the purely combinatorial problem of 
finding the group factors corresponding to the Wick contractions. What is left is 






where c 2k (N) is the sum over all possible traces of 2k T a matrices suitably contracted. 
Fortunately this group factor is generated by a matrix integral [2T|, and the final result 
reads 
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Wn[A] 



N 



exp 



g 2 n z A 



■W ^2 2 



Lj^foV.A/2) 



(38) 



reproducing the zero-instanton result Eq. ( |23"D obtained via non-perturbative methods. 

B. Resummation of the perturbative series defined via 't Hooft prescription 

As already hinted, we were not able to resum the full perturbative expansion for W n pre- 
scribed with 't Hooft. We believe after exertion that such a formidable task is, if not 
impossible, at least extremely tough. Nevertheless we were successful in the computation of 
the Wilson loop with n windings at 0(g 4 ), which suffices to give a flavour on how things work 
at higher orders. Moreover we performed the limit n — > oo, n 2 A fixed, and, by exploiting 
the duality n *-* N, the large-N limit, with g 2 N fixed. 

Firstly, let us start from the perturbative definition of W n in the light-cone gauge (A- = 
0) 

W n = -^AfTr | fvA^ exp i f d 2 x (--F^F» va + J^A 111 

5{A a _) V exp (ig&Al T a dxA 1 , (39) 

evaluated at J a = 0. It easy to recognize that Eq. ( |39"D can be rewritten as follows 



W n =^AfTr { Vvxp 



9<pT a 



5J a (x 



■ dx~' 



exp 



\ d 2 xd 2 y.J a (x)D(x-y)J a (y) 



j=o 

(40) 



where now the propagator D(x — y) is defined through the CPV prescription Eq. (^). We 
consider a light-like rectangle with sides 2L, IT (see Fig. 1) and choose the currents with 
support on the contour, so that 

J a (x + ,x-)=j a u 6(x--L) +j a d S(x-+L). 



With this choice the perturbative expansion Eq. ((5]) for W n reads 



W n = ^AfTt{Vvxv> 



■V exp 



Lr f ^ 



dx~ 



\ xr> 



exp 



T 



dx A 



V exp 
Pexp 






£ 



T-IO 



*#(* + ) 

<5 



dx 



Ci, 



W* + ) 



ete 



iL 



dx + ^(x+)^(x + ^ 



(41) 



3=0 



with C*i and C*2 as in Fig. 1. Clearly in the limit j£, j% — > the only non- vanishing 
contributions are those with a matching number of derivatives with respect to j u and j^. 

Now everything is settled and we can easily derive the expression of Wn , i-e. the Wilson 
loop with n windings 0(g A ). The following prefactor is common to all classes of diagrams 

9 4 A 2 



g*(iL) 2 {2T)< 



(42) 



A = 4LT being the area of the loop. As a matter of fact, g 4 comes from four derivatives 
with respect to the currents which contribute at this order, (iL) 2 is produced when two 
derivatives act on the last exponential in Eq. (f41~l) , which represents the only non- vanishing 
contribution, and finally (2T) 2 is given by integration over the loop variables. Notice the 
expected dependence on the area due to invariance with respect to area-preserving diffeo- 
morphisms typical of d = 2. 

In Appendix A we show which classes of non-crossed and crossed diagrams contribute to 
Wn and also provide the exact counting (as a function of n) and area factor. The latter 
turns out to be either \ or 1 depending on the presence of integrals in the loop variables 
which are nested as a consequence of the definition of the V exponential. Furthermore one 
has to take into account that non-crossed diagrams produce an additional factor N 3 /4 from 
Tr \T a T a T b T 6 ] , whereas the analog for crossed diagrams is N/4 from Tr \T a T b T a T 6 ] . 
Remember eventually the factor 1/N appearing in Eq. ([4l"l). 

We present here the final results, which read 

9 4 A 2 n\ T2 






96 

g 4 A 2 n 2 

96 



A^ 2 (2n 2 + 1) 
(n 2 - 1) 



(43) 
(44) 
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and, summing up, 



W« = - 9 -^r~ ( 2n2N2 + iV 2 + n 2 - 1) , (45) 

9o 



which coincides with the coefficient 0(g 4 ) of the expansion of Eq. ([H|) obtained via non- 
perturbative methods. Moreover, already at this level, the symmetry under the exchange 
n <-> N is manifest once the area is rescaled by N/n. 

Let us now go back to Eq. flil) ) in order to extrapolate the relevant limits n — > oo and 
N — > oo. We have already stressed that they are dual to each other (provided the loop area 
is suitably rescaled), so that we will explicitly perform only the former and deduce the latter 
by symmetry. In Sect. II the large-n limit of the non-perturbative W„, with An 2 fixed, was 
proved to coincide with the resummation of the perturbative series when WML prescription 
was adopted and A — ► An 2 . Thus we guess that at large n dominant configurations of 
loop diagrams where 't Hooft prescription is imposed must be weighted by the same area 
factor. This precisely happens for WML as we realized in Sect. Ill A. In fact, by inspecting 
the computation of W n at order 0(g 4 ) we are easily convinced that in this limit dominant 
configurations, either crossed or non-crossed, are those which spread through the maximum 
number of sheets of the n-fold loop, i.e. those in which the vector fields are all attached 
on different sides of the loop. For instance the leading terms at 0(g A ) are those shown in 
Figs, li, 2d. Going further, at a generic order g 2k one needs at least Ik sheets in order to set 
up such configurations, and their number is straightforwardly determined. Actually, there 
are ( 2 ^,) ways of extracting 2k sheets out of n and in addition we have to take into account all 
possible configurations of the basic 2fc-sheet structure with any kind of cross-over. Therefore, 
the counting of dominant configurations is given by 

TV 

x c 2k (N) , (46) 



(2k)\{n-2k) 

where C2k{N) is the coefficient we acquainted with in Sect. IIIA. Hence, at the leading order 
in n, we end up with 



«4£(-^)^ 



2 \ & 

\ 2 

fc=0 
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and, switching to the euclidean space via a Wick rotation 



Wn-HxiA) = T7 ex V 



4^ ( ^ ) , (48) 



-\9 2 An 2 

4 

which corresponds to Eq. (p3|). 

The next step consists in performing the large- N limit of the resummed perturbative 
series in the axial gauge with the light-cone CPV prescription. The duality n <-> N allows 
us to conclude 



Wiv^ooOA) = -exp 
n 



1 ~2 A 

--g 2 An 
4 



i& ( '-^ ) , 



where the rescaling in Eq. fl26|) is understood. Apart from the numerical coincidence, the 
duality entails a deep relation among different faces of the same theory. Naively one could 
think that the large- N limit, being dominated by non-crossed graphs in which the trace 
factor considerably simplifies, is easy to recover . Nevertheless, a surprising complexity in 
the integrals over loop variables arises, since every class of diagrams is characterized by a 
different area factor. This is already manifest from the computation of the Wilson loop 
at order g 4 (see Appendix A). Thus the most relevant conclusion we can draw is that an 
enlightening interplay between geometry and algebraic invariants takes place, connecting 
the large-TV with the large-n limit and, in a sense, the strong and the weak phase of the 
theory. 

IV. THE SPECTRAL DENSITY 

U(N) Yang-Mills theory on the sphere S 2 exhibits in the large- A limit with g 2 N fixed, a 
third order phase transition (Douglas-Kazakov (DK) phase transition) [p4f] . According to 
the total area A of the sphere, two different regimes occur: a strong coupling regime when 
A > ix 2 (in suitable units, i.e. ^— = 1) and a weak coupling regime with A < n 2 . The 
difference between them can be traced back to the presence of instantons on the sphere 
23| , |2B|P^] , which fully contribute for large A, but are instead exponentially suppressed in 



the weak-coupling phase. 
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The two regimes can be conveniently characterized by the density p of Young tableaux 
indices which single out the various irreducible representations of the group U(N). In 
the strong coupling regime the distribution p undergoes a phenomenon of "saturation" in 
the large-iV limit, namely it is equal to one in a finite interval of its domain. In turn, if 
we consider a Wilson loop on the sphere with the contour C running around an equator 
and winding an arbitrary number n of times, the density p is related by a kind of duality 
transformation |22| to the spectral density a on the unit circle of the eigenvalues of the 



unitary operator Uc, representing the loop. 

The strong-coupling phase is characterized by the fact that those eigenvalues completely 
fill the unit circle, whereas in the weak coupling regime a gap develops where the distribution 
a vanishes. Before decompactifying the sphere to the plane by taking the limit A — > oo, we 
have to learn how the density a evolves when we replace the equatorial contour with a smooth 
arbitrary one. In ref. p2] it is shown that, when one of the two different areas generated in 
this way, let us call it A, becomes small enough, a gap in the eigenvalue distribution appears 
even in the strong DK phase (i.e. for A > it 2 ). Actually, when the total area A approaches 
its critical value it 2 from above, A becomes critical at its maximum value A cr = ^-, namely 
when the contour is along an equator. 

When the sphere is decompactified to the plane by taking the limit A — ► oo at fixed A, 
we are automatically in the DK strong phase; it is then remarkable that A cr reaches, in the 
limit A — > oo, a well defined value, A cr = 4. 

We are going to obtain all these results from the explicit expressions we got in Sect. II for 
Wilson loops with an arbitrary winding number W n . We shall explicitly discuss the large- N 
case on the plane (A — > oo). The spectral density cr(A, 9), A being the area enclosed by the 
loop, is related to W n by the following Fourier series |Tj 



a(A, 9) = i 

Z7T 



l + 2^cos(n0)W n 



(50) 



n=l 

We shall compute this function and show that, according to different values of A, a gap 
may develop in the ^-variable in which the distribution vanishes. This phenomenon was 
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first noticed by Durhuus and Olesen |I|, who derived the spectral function solving a kind 
of Makeenko-Migdal equation with suitable boundary conditions. 
By introducing Eq. fl3~T|) into Eq. floTf) we get 



<?(A,9) = ±- 

Ztt 



oo ^ 



n=l 



where 



F(A 



oo _. 



n 



-n(4-i0) 



^-[l + 2ReF(A,9)], (51) 

Z7T 



n=l 



A e -^fi + i 



(52) 



* r 27ri y t 

Here we have used a well-known integral representation for the Laguerre polynomials in 
which the integration contour T encircles the origin of the complex t-plane. 
When the inequality 



-Mt+\) (l + -\ 



< 1 



(53) 



is satisfied, the series in Eq. (|52"D can be summed giving rise to a single-valued function 



y- 1 e - n [A(t+i)-ie] ( l + l 
^ n \ t 



n=l 



log 



-A(t+^)+i6 



1 

1 + 7 



(54) 



The logarithm in turn can be analytically continued; as is well known, it has branch points 
when its argument either vanishes or diverges. Whereas it is easy to realize that the only 
divergence occurs, at finite \t\, at t = 0, a careful study of the roots of the trascendental 
equation 



is in order before drawing any conclusion. 

It is amusing to notice that the change of variable 

exp(t£) 
1 - exp(z£) 

turns Eq. (0) into Eq.(3.4) of ref. |B|. 

To the success of our calculation it is essential that the contour T in Eq. 
deformed while the logarithm still remaining single-valued on it. 



(55) 



(56) 



can be 
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The discussion of the roots of Eq. ( |55| ) is deferred to Appendix B. There we show that, 
for A > 4, only two roots are possible and only one of them is encircled by the contour I\ 
If we set t = x + iy and denote this root by the values (x, y), which obviously are functions 
of A and of 9, a cut can be drawn along the segment (0,0) — {x,y) in the complex t-plane 
and the contour T can be deformed to the contour 7, encircling the branch points and just 
running below and above the cut. It is then immediate to conclude that 



F(A, 9) = x + iy 



and that, eventually, 



*{A,9) = ~(x + ±). ir,7) 



As long as A > 4, x + | is positive for any value of 9. The eigenvalues are filling the 
entire interval (— 7T, it). 

When A < 4 the situation dramatically changes. Although we are in the strong DK 
phase (we are on the plane), a gap starts showing up in the eigenvalue density for 9 close to 
7r ( and to — 7r; remember a is an even function of 9). Actually, for A < 4, there is a value 
9 cr < tc such that x(A, 9 cr ) = — ~. Once this value of 9 is reached, x remains equal to — |, 
and consequently a vanishes. 

In Appendix B we show that the expression for 9 cr coincides with the one in Eq.(4.8) of 
ref. |16|, namely 



9 cr = \IA h arccos ( 1 ) . (58) 

Eventually, in the limit A — ► 0, the density approaches the periodic (^-distribution, as it 
should. 

Actually we have already seen in Sect. II that, in the small- A limit, instantons on the 
sphere are suppressed even in the strong DK phase. This phenomenon has a simple geomet- 
rical interpretation: it occurs when one of the two areas singled out by the loop gets small 
when compared to the instanton size. It might then seem that the different behaviour of 
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a according to the values of A, is fully driven and characterized by instantons. This issue 
would in turn be of crucial relevance in singling out peculiar properties of the light-front vac- 
uum when compared to the one in equal-time canonical Fock space after decompactification 
of the sphere to the plane. 

Unfortunately the situation is not as simple as that. In ref. [Q the contribution to the 
Wilson loop from the zero-instanton sector on S 2 was clearly pointed out. It was also noticed 
that, in the decompactification limit, this contribution exactly coincides with the sum of the 
perturbative series when graphs are evaluated according to the WML prescription for the 
propagator. In Sect. II we have generalized this result to a loop with an arbitrary winding 
number n. It is then a fairly simple exercise to take the large-iV limit and to insert the 
result in Eq. (|50|) 

"" l 1 + 2 

n=l 

J\ being the usual Bessel function. Of course its geometrical meaning is now completely 
different. 

Still, due to the absence of instantons, one would naively believe that a ML should behave 
like in a weak phase, in particular should develop a gap. But this is not always the case; as 
a matter of fact, using a standard integral representation for the Bessel function, we get 

a ML = - [ dtVT^5 P (9 + 2tVA), (60) 

5p being the periodic 5-distribution, leading to 



Ecos(n9) T . /—. 



-7T < e < 71- , (59) 



mt 1 \-v / (2-nm + e) 2 n , . 

the sum running over all positive and negative values of m such that the argument of the 
square root is positive. 

It is easy to check that this expression reproduces the limits 

lim a ML = — 

.4^00 2ir 
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and 



\ima ML = 5(9). 



For A < ^j- there is indeed a gap, similar to the one in the weak phase on the sphere, 
but with a new threshold. For larger values the gap disappears. 

Only the functional form (a square root) is reminiscent of a weak phase. There are 
moreover discontinuities in the derivative of the density anytime a new value of m starts 
contributing. 

One would eventually conclude that the occurrence of a gap is not directly related to the 
presence of instantons on S 2 ; even in the zero instanton sector, after decompactification, no 
gap occurs for sufficiently large values of A. Only when A is small enough, the complete and 
the zero-instanton solution exhibit similar behaviours. We emphasize that this statement 
has to be proven, since the limits N — > oo and A —>■ could in principle not commute. For 
small area, only the term with m = contributes to the summation Eq. (|6lD , and therefore 
we deduce 9^ L = 2y/A. This is precisely the leading term in the expansion of Eq. fl58|) , so 
that the two values of 9 cr coincides at small area. Let us now show that also a (A, 9) Eq. (|57|) 
and a ML (A, 9) have a similar behaviour for small A. As 9 < 9 cr , by defining 9 = X9 cr , we 
can write 



a 



ML 



VT^X 1 . (62) 



nVA 

Notice that the last factor in Eq. (B3) is finite as A — > 0. By comparing Eq. ( |62| ) and Eq. (p7|). 
we expect that the solutions (x, y) of Eq. fl55|) behave for small A as {—\ + -4=, -jj), a- and b 
being finite numbers. As a matter of fact, solving the trascendental equation Eq. ( |o^| ) with 
this ansatz in mind, we find 

V - ~Tt ■ ( 63 ) 



Thus, we come to the remarkable result 
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r. = - ( x+ + I ) ~ -^VT^X 1 , (64) 



which is just Eq. fl6~2f) . 

Its precise meaning is elucidated if we remember that cr(^4., 0) has a distribution character 
in the limit .4. — >• 0. After introducing the test function <p(9), and defining 

/9cr 
d9ip(9)a(A,9), (65) 

■u cr 

we get 



.4^0 



(a,<p)~- d\VT^X*<p(\6 cr )^<p(0), (66) 

which shows indeed that, when .4. — » 0, cr(A, 9) —* 5(9) in the topology of distributions, 
both for the exact and the WML case. 

On the other hand, at large areas, the genuinely perturbative solution (i.e. the one 
coming from the equal-time quantization and corresponding to the zero-instanton sector), 
becomes less and less reliable. 

To gain a better insight, we consider again the zero-instanton sector on S 2 , before decom- 
pactification. Eq. (|2^) teaches us that this is easily achieved by replacing A with A ■ 
In the equatorial condition A — A = A, Eq. (|6T|) becomes 



^^Wi-e^, -»^«- «") 



ny/A^ V A 

It is amusing to notice that the critical value for A is now again A cr = ti 2 and that, below 
such a critical value, the ML spectral distribution exactly coincides with the one in the weak 
DK phase on the sphere p4}] . 

Actually, from Eq. (|22|), we can derive a stronger result; we can indeed follow the evo- 
lution of the critical area of the loop as a function of the size of the sphere also for the 
zero-instanton sector. As a matter of fact, Eq. ( |6~H ) becomes 



a ML (A^ L ,tt) = - VA J I - AAMJ f A A AMT „ (68) 

ny/A% L (A-A% L )y 4A% L (A-A% L Y 
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n 2 

A 



(69) 



leading to 

A 
aML __ 2_ 

cr ' 2 

Starting from the value A^. L = \, when A = 7r 2 , the same as the one of the exact solution, 
in the decompactification limit A — > oo we recover the threshold A^ L = 7, lower than the 
one of the exact case (A cr = 4). 

V. CONCLUSIONS 

The main concern of this paper was to gain a better insight in the dynamics of YM 2 by 
considering the class of Wilson loops winding n times around a closed smooth contour, 
either on a compact manifold (S 2 ) or on the plane. As already noticed by several authors, 
thanks to the invariance under area-preserving diffeomorphisms, these loops carry enough 
information to characterize the theory in the large- N limit. 

Our deepest result is probably Eq. (|26|). It shows a curious interplay between two integral 
numbers, iV characterizing the internal symmetry group and n representing the number of 
windings along the contour. In a geometrical language, the latter controls the displacement 
of the connection on the base manifold, the former is related to displacements along the group 
fiber. In two dimensions, the winding n has essentially a topological character, as it appears 
in the n-th power of the group element Uc\ the Wilson loop can then be interpreted as 
a generalized Clebsch-Gordan coefficient intertwining different U(N) representations. This 
is, we believe, the root of Eq. ( ^6|) and, of course, it entails far-reaching consequences. In 
particular it relates the abelian-like exponentiation for U(N) when n = 1 to the genuine 
triviality of the Maxwell theory (N = 1) for any n. 

All these features are reproduced in a perturbative context. Actually concrete exam- 
ples of different colour and geometrical interplays are explicitly exhibited, according to the 
different expressions used for the vector propagator (either WML or 't Hooft). When the 
't Hooft propagator is used the perturbative series leads to the exact result. However, for 
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n > 1, the Wilson loop feels the non-abelian nature of the theory. The winding number n 
probes its colour content. The related light-front vacuum, although simpler than the one in 
the equal-time quantization, as it automatically takes instanton contributions into account, 
cannot be considered trivial. 

A more intriguing result concerns the limit n — > oo keeping n 2 A fixed, namely the 
small-area limit. In such a limit (Eq. (|3^)), one exactly recovers the contribution of the 
zero-instanton sector, namely Eq. fl23l) . Intuitively this can be understood by remember- 
ing that instantons possess a typical length scale measured by their size. One might then 
be tempted to conclude that the phase transition at large N is completely driven by in- 
stantons. Unfortunately this characterization is true only for small areas, when instantons 
cannot contribute. In Sect. IV it is shown that, for areas larger than a critical threshold, the 
spectral function which characterizes the eigenvalue distribution in the exact solution, when 
computed retaining only the zero-instanton sector, does not exhibit any gap. At large areas, 
the genuinely perturbative solution (i.e. the one coming from the equal-time quantization 
and corresponding to the zero-instanton sector), becomes less and less reliable, as expected 
on a general ground. It is then quite remarkable that the perturbative light-front resum- 
mation succeeds in reproducing not only the correct string tension, but also the non-trivial 
corrections (see Eq. (|17])) due to colour- winding factors, the winding number n carrying the 
information about the non-abelian topology. 

One should perhaps conclude with a comment concerning higher dimensional cases, in 
particular the dimension d = 4. We believe that most of our results are typical of the 
two-dimensional case. Perturbation theory at d = 2 + e is discontinuous in the limit e — > 
|5|j27[| ; on the other hand the invariance under area-preserving diffeomorphisms is lost when 



d > 2. In a perturbative picture the presence of massless "transverse" degrees of freedom 
(the "gluons") forces a causal behaviour upon the relevant Green functions, whereas in the 
soft (IR) limit they get mixed with the vacuum. The light-front vacuum, which also in two 
dimensions is far from being trivial, in higher dimensions is likely to be simpler only as far 
as topological degrees of freedom are concerned. Of course there is no reason why it should 
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coincide with the physical vacuum since, after confinement, the spectrum is likely to contain 
only massive excitations. Moreover, to be realistic, "matter" should be introduced, both 
in the fundamental and in the adjoint representation. Therefore, before going to higher 
dimensions, our two-dimensional considerations should perhaps be generalized to the case 
in which "matter" is present. Although many papers have appeared to this regard in the 
recent literature, we feel that further work is still needed to reach a complete understanding. 

VI. APPENDIX A: W n 0(g 4 ) WITH 'T HOOFT PRESCRIPTION 

In this appendix we explicitly compute W n at order 0(g A ). We have seen in Sect. Ill 
that at this order all diagrams come with the prefactor given by Eq. (0). We now have 
to group all possible diagrams in different classes according to their contribution either to 
Tr \T a T a T b T fe ] (non-crossed diagrams) or to Tr [ T a T b T a T 6 ] (crossed diagrams), and to 
the area factor. We then have to specify how many of them belong to each class (depending 
on the number of windings n) and the area factor they produce []. Such numbers are 
summarized in the following tables: Table [I] refers to non-crossed diagrams depicted in Fig. 2, 
whereas Table || to crossed ones depicted in Fig. 3. Finally, summing up the contributions 
from all loops, each multiplied by the prefactor (|42|) , the proper trace and area factor, and 
taking into account the multiplicity, we recover the result announced in fl45|) for the Wilson 
loop with n windings at order 0{g i ). 

VII. APPENDIX B 

In this appendix we will sketchily show that Eq. fl55|) admits two solutions and, of these, 
only one can be encircled by the integration contour. Let us start from Eq. (|55|), which, 



4 Recall it is either ^ or 1 depending on whether there are nested integrals in the loop variables 
or not. 
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after inserting t = x + iy, reads 

1 + X + l l e -A (x+l) e i{6-Ay) = 1 _ / 70 x 

x + iy 
The first obvious symmetry to be noticed is (x,y,9) —>■ (x, —y, —9). In addition, by taking 
the absolute value of Eq. (|70|) , the symmetry x — > —x — 1 is also manifest. Combining 
Eq. (|70| ) and its complex conjugate, we obtain an equation for the shape of the boundary 
which separates the region of convergence of Eq. ([52]) and the forbidden region 

y 2 = -(x 2 + x + \) + (x + \)cot\i{A{x + \)), for x±-\. (71) 

It follows 

~te = "~ 'sinh 2 (A(x+±)) fiX ' ,A: '' (?2) 

where 

fix; A)=A- SinM 2 f^ 1 +1)) + 2 sinh 2 (A(x + ±)) . (73) 

It is easy to check that the same function / appears in the derivatives of the coordinates x 
and y with respect to 9. From Eq. ( [70|) we infer 

dx 8yyg(x,y) 



89 89 f{x; A) ' 
8y ffaA) 



86 f 2 (x;A)+y 2 g 2 (x,y) 

with 

1 

9 fay) 



(74) 



(x + l) 2 + y 2 x 2 + y 2 
By comparing Eqs. (|72|) , (|74]) we deduce 



dy cfa/ 2 



sign — = —sign -j— it x > — I , 



sign M = sign ^ if x< - ±2 - (75) 

We conclude that, if for 9 = 9 a root of Eq. (|S5|) exists, for 9 > 9 its ordinate ?/,-(#, ^4.) 
follows the boundary continuously. Moreover, Eq. (^) tells us that for 9 = there are two 
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roots sitting on the x-axis (y r = 0), whose x coordinates are the solutions of the following 
equation 

x + 1 A(x+V 



2, 



X 



(76) 



One can check that Eq. ( |7E|) admits only two solutions, for any value of A, symmetric with 
respect to the axis x = — \. Clearly, when x approaches the value — |, the condition y 2 > 
in Eq. (|n]) can be fulfilled only if A < 4. Therefore, we start at 9 = with two roots on 
the real axis which move continuously along the boundary as 9 increases if A > 4, whereas 
for A < 4 there is a critical value of 9, let us call it 9 cr , for which the two roots reach the 
axis x — — \ and then remain on it, moving in the imaginary direction. Hence a gap is 
originated, since 

a(A, 9) = for 9 > 9 cr , A < 4 . 

At this stage we present an hand-waving argument to show that only one of the two roots 
is encircled by the integration contour. Nevertheless a rigorous proof based on a thorough 
examination of Eq. (p3|) can be given. Although the region of convergence of Eq. (|52D 
appreciably varies according to the value of A, it is always made up of disconnected pieces 
and the axis x = —\ represents a border line between a domain of convergence and a 
forbidden region. As a consequence, the integration contour, which encircles the origin of 
the complex plane (lying in the forbidden region), is forced not to cross such an axis and 
to capture just one of the two roots (remember they are symmetric with respect to the 
exchange x — > —x — 1). This root in turn and the origin become the two branch points of 
the logarithm in Eq. (|54]). 

Let us now determine the exact value of 9 cr in case A < 4. The criticality is reached 
when x = — \, so that we evince from Eq. (|7T|) 



-H 



Then Eq. (|70| ) becomes 
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1 - l cos (6> cr - Ay cr ) - y cr sin (6> cr - Ay, 



2-2 v - —, ^ v - „-, (77) 

y cr = i sin (9 cr - Ay cr ) + y cr cos (9 cr - Ay cr ) 
Eventually, by inserting the value of y cr , the solution is straightforwardly found 



e cr . = \jA h arccos ( 1 ) . (78) 
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TABLES 



Diagrams 


number 


area factor 


Fig. la 


n 


l 

2 


Fig. lb 


2n(n - 1) 


1 
2 


Fig. lc 


n(n — 1) 


1 
2 


Fig. Id 


n(n— 1) 
2 


1 


Fig. le 


n(n— 1) 
2 


1 


Fig. If 


n(n — l)(n — 2) 


1 
2 


Fig. lg 


n(n — l)(n — 2) 


1 


Fig. lh 


n(n— l)(n— 2) 
2 


1 


Fig. li 


n(n-l)(n-2)(n-3) 
3 


1 



TABLE I. Classes of non-crossed diagrams contributing to W n 0(g 4 ) with counting and area 
factor. 



Diagrams 


number 


area factor 


Fig. 2a 


2n(n - 1) 


l 

2 


Fig. 2b 


n{n — l)(n — 2) 


1 
2 


Fig. 2c 


n(n— l)(n— 2) 
2 


1 


Fig. 2d 


n(n-l)(ra-2)(n-3) 
6 


1 



TABLE II. Classes of crossed diagrams contributing to W„ 0(g ) with counting and area factor. 
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FIGURES 



FIG. 1. Contour. 




(a) 




(b) 



35 




(c) 





(d) 




(e) 






(g) 
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(h) 



FIG. 2. Non-crossed graphs. The i th sheet corresponds to the i th winding. 





(b) 
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(c) 




FIG. 3. Crossed graphs. 
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